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Based on Takahashi-Umezawa thermo field dynamics and the order-invariance of Weyl ordered 
operators under similar transformations, we present a new approach to deriving the exact Wigner 
functions of thermo number state, photon subtracted and added thermo vacuum state. We find that 
these Wigner functions are related to the Gaussian-Laguerre type functions of temperature, whose 
statistical properties are then analysed. 



I. INTRODUCTION 

In recent years photon subtracted and added quantum states have been paid much attention 
because these fields exhibit an abundant of nonclassical properties and may give access to a complete 
engineering of quantum states and to fundamental quantum phenomena [1-8]. However, all these 
discussions are restricted to the case at zero point temperature. In fact, most systems are not 
isolated, but are immersed in a "thermal reservoir" , excitation and de-excitation processes of a 
system are influenced by its energy exchange with reservoirs. In this work we study field properties 
by photon subtracting and adding at finite temperature. 

The Wigner function (WF) is a powerful tool to investigate the nonclassicality of optical fields 
[9,10]. Its partial negativity implies the highly nonclassical properties of quantum states and is often 
used to describe the decoherence of quantum states [7,8,11,12]. In one dimensional case, the WF of 
a density matrix p is defined as Tr [pA(a)] , where A.{a) is the single-mode Wigner operator, whose 
normally ordered form and Weyl ordered form are given as [13-15], respectively, 

A (a) = - : e-(?-Q)'-(P-^)' : ^ i : e-2("-)("*-^) : , (l) 

TT TT 

and 

A{a)^^''S{a-a)d{a* ~a^)[, (2) 

where a — {q + ip) /\/2, a = {Q + ±P) /V^, [Q,P] = i, h — 1; a and a'^ ([a, a^] — 1) are Bose 
annihilation and creation operators, the symbols : : and ; : denote the normal ordering and the 
Weyl ordering, respectively. Our main aim is to provide a new and direct approach to deriving 
the WFs of quantum states at finite temperature by using the order-invariance of Weyl ordered 
operators under similar transformations [13-15], which means 

S] (ooo) 'S-i = ■S'(ooo)S'-i;, (3) 

as if the "fence" [ ', did not exist, so S can pass through it. We also appeal to the Takahashi-Umezawa 
thermo field dynamics (TFD) [16-18], we consider it convenient to obtaining the explicit expressions 
of WFs. 



II. BRIEF REVIEW OF THERMO STATE 



The main point of TFD lies in converting the evaluation of ensemble average at nonzero tem- 
perature into the equivalent expectation value with a pure state. This worthwhile convenience is at 
the expense of introducing a fictitious field (or a so-called tilde-conjugate field, denoted as operator 
5^) in the extending Hilbert space H, thus the original optical field state \n) in the Hilbert space 
Ti is accompanied by a tilde state \n) in H. A similar rule holds for operators: every annihilation 
operator a acting on Ti, has an image a acting on H. At finite temperature T the thermal vacuum 
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|0(/3)) is defined by the requirement that the vacuum expectation value agrees with the statistical 
average [16-18], i.e. 

{A) = Tr ip^A) = {OmA |0(/3)) = Tr {Ae-^'^) /Tr {e-^") , (4) 

where f3 = -jyp, k is the Boltzmann constant and H is the system's Hamiltonian. For the ensemble 
of free bosons with Hamiltonian Hq — uja'^a, the thermal vacuum state |0(/3)) is 

|0(/3)) = sech^exp [a^5^tanh6i] |0, O) = S" (61) |0, O) , (5) 

where |0,0) is annihilated by a and a, [a, a^] — 1, and 

S (e) = exp [61 (fltfl^ - a~a)\ , (6) 

is the thcrmo squeezing operator which transforms the zero-temperature vacuum |0, O) into the 
thermo vacuum state |0(/?)) , and 6 is related to the Bose distribution by 

tanh^ = exp(-^), (7) 

which is determined by comparing the Bose-Einstein distribution 

(8) 
and 

(0(/3)| a^a |0(/3)) ==sinh2 61. (9) 

In particular, when operator A is the Wigner operator A [a) itself, it is easy to see that 

Tr4A(a)e-'^^)/Tr„(e-^«) - (0(/3)| A (a) |0(/3)) 

= Tr„,a[A(a)|0(/3))(0(/3)|], (10) 

which is just the WF of thermo vacuum state. From Eg. pTI)) one can see that the calculation of 
WF for thermo states is converted into the expectation value of Wigner operator in themo vacuum 
state |0(/3)) (pc |0(/?)) (0(/?)l)i which is defined in the enlarged Fock space. This implies that it 
is convenient to deriving some WFs of density operators at finite temperature by doubly enlarging 
the original space. 



exp 



(^ 



III. NORMALLY ORDERED FORM OF S'^ (6) A (a) S (6) 

In order to deriving conveniently the WFs of density operators at finite temperature, let's first 
calculate the normally ordered form of (6) A (a) S (9) . Recalling that for single-mode case the 
Weyl rule [13-15] is defined as 

H{a,a^) =2 j d^ah{a,a*)A{a), (11) 

where h {a, a*) is the classical function corresponding to operator H (a, at) . Eq. pT]) expresses the 
Weyl correspondence rule, using ^ it can be expressed as 

H[a,a') = J d'^ah{a,a*) S {a — a) 6 [a* — a^) 

= [h{a,a^) (12) 

which means that Weyl ordered of operator ; h (a, a^) ; , whose Weyl correspondence is h {a, a*), can 
be obtained by just respectively replacing a, a* in h {a, a*) by a and without disturbing the form 
of function h. 
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According to the Weyl ordering invariance under similar transformations [13] and the following 
transform relation 



it is easily seen 



S^e)aS{9) = acoshe' + atsinh6l, 
S''{9)dS{e) = acosh6' + atsinh6i, 



S'' {0)A{a)S{e) = - S{a~acosh9-a''smhe) 
x5 (a* — cosh 6* — asinhS) 



(13) 



(14) 



which is just the Weyl ordering of {6} A (a) S {&) in the enlarged Fock space. Based on the Weyl 
rule, the classical correspondence h (/?, p* ; /3, /3*) of the operator S^t (6) A (a) S [9) can be obtained 
by replacing (a, a^) and (a, fit) with {(3,l3*) and {(3,(3*), respectively, i.e., 

h(^(3,(3*;(3,(3*^ = i(5 (a - /3cosh6' - /3* sinh6') 

(15) 



(16) 



x6[a* -P* cosh6'-/3sinh6'] 
It then follows from Eqs. (fTT|) and p5|) that 



where A (^(3 , (3* ; (3 , (3 

A(/3,r;/3,/3*) 



is the two-mode Wigner operator, whose normally ordering form is 
1 



exp 



(at„^*)(a_^)_2(at-/3*) (fi-/?) 
On substituting Eq. (fT7| into Eg. p^ and using the integral formula [19] 

/ ^e^l^l'+«-+''^* = --e-f, Re(C) < 0, 
7 71- C 

we can derive the normally ordered form of (|16p as follows 

S'' i9)A{a)S{9) = 2 y" ^^^^(5 (a - /3cosh6i - /3* sinh^) 

xS(a* - (3* cosh 9 - (3 sinh 9^ 

(flt - p*) {a-P)-2 (fit - (3*^ (fi - /?) 



(17) 



(18) 



X : exp 
sech26' 



-2|Qpscch2e . 



exp |— (afi + fltfit) tanh 29 

-|-2sech26' [sinh 6* (Q;*fit -(- afi) + cosh 6* {a* a + aa)) 

- (fitfisinh^e + fltacosh^e)]}: , (19) 

which is just the normally ordered form of (|16p . Eq. (fT9|) directly leads to the WF of the thermo 
vacuum state |0 (/?)), 



(0(/3)|A(a)|0(/3)) = (O,6|5t(0)A(a)5(^?)|O,O) 



sech26' 



TT 



-2|Qpscch2e 



1 - e"''" 
7r(TTe-^' 



-2\ay 



(20) 



IV. WIGNER FUNCTION OF PHOTON-SUBTRACTED THERMO VACUUM STATE 

At finite temperature, the photon-subtracted thermo vacuum state can be expressed as [20] 

Pi = Cia"|0(/3))(0(/3)|at", (21) 
where Ci is the normahzed factor, defined by 

Cf 1 = Tr [a''S (6) |0, O) (O, 0| (9) a^"] , (22) 
which can be calculated as follows. Using Eq.© and the binomial formula 



1=0 



we have 



Cf^ = (0, 6| 5^ (6*) at"a"5 (6*) |0,0) 

= sech20(O,O|e''"*''"''^a^"a"e''''^^*^"''^ |0,6) 



sech^ei J2 tanh'+'=6i/A:,fc 



k.l=0 



IJ 



= sech^g V „ • , tanh^' 9 = nl sinh^" 0. (24) 

By using Eqs. (PT|) and (|19p . we calculate the WF of photon-subtracted thermal state pi 

Wiia) = Ci (0,0| 5^ (61) a^"A (a) a"S' (61) |0,0) 

= (0, 0| [S^ (9) a^"^ (9)] (9) A (a) S (9) [S^ (9) a"5 (9)] |0, 0) . (25) 



Noticing Eq. ()T3|) we see 

[S'' {9)a''S{9)]\0,b) = (acosh6i + a^sinh6i)" |0,0) 

= \/^sinh"6i|0,n) , (26) 

then substituting (|26p into Eq. ((25|l and using Eq.(fT9|) yields 

Wi(a) - ^ (nle^" (sech20)" "e^""|n) 

TT cosh 26* 



g-2|a|^soch2e « ^*k^i /2sinh6l 



TT cosh 26* ^ fc!;! V cosh 26* 



k+l 



^ b;r7^ (n|at'=(sech2^)-aM") 



n cosh"+^ 29 f^^nnin-iy. V cosh 20 



7171 rr, _ 1\\ rn^hOf) / ' ^ ' 



where we have used the identity operator [21] 

exp [Aa^a] =: exp [(e^ - l) d^'a] : . (28) 
Recalling that the definition of Laguerre polynomials [22] , 

L„{x)^y {-xY, (29) 

toil'-) i^-iy- 



4-4 
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FIG. 1: Wigner function distributions of photon-subtracted thermo state in {q,p) phase space with (a) 
n = l,e = 0.2, (h) n = l,e = 0.8, (c) n = 2, = 0.8, and in (|a| , 6) space with (d) n = 1. 

Eq. (j27p can be further put into the following neat form, 



which is just the WF of photon-subtracted thermo vacuum state, a Gaussian-Laguerre type function 
of temperature, since tanh^ = exp (— 2^^)- Due to cosh 20 > and 7^17^ lo^l^) ^ 0, for the 

photon-subtracted case, Wi (a) has no chance to present the negative value in phase space, which 
can be seen from Fig.l. On the other hand, the amplitude value of WF in {\a\ , 9) space decreases 
with the increasing temperature (corresponding to 6). In appendix A, in order to check the result in 
Eq. (I30|) . we have derived the WF of photon-subtracted thermo vacuum state by using the coherent 
state representation of Wigner operator. Comparing with the result in Ref.[20], Eq. (|30p seems more 
concise and convenient for further discussion. 



V. WIGNER FUNCTION OF PHOTON-ADDED THERMO VACUUM STATE 

At finite temperature, the photon-added thermo vacuum state is expressed as [23] 




(30) 



P2 = C2a^"|0(/3))(0(/?)| 



a 



.n 



(31) 



By the same procedures as deriving Eqs. (l22|l and ((26)) . we have 

C2~^ = n! cosh^" 6*, 



(32) 



and 



S'^ {9) a1'"S' {9) |0, 0) = %/^cosh" 9 \n, O) . 



(33) 



FIG. 2: Wigner function distributions of photon-added thermo state in {q,p) phase space with 9 — 0.2 for 
(a) n = 1, (b) n = 2, and in (|a| , 9) space with (c) n = 1 and (d) n = 5. 



a Gaussian-Laguerre type function which may present negative region in phase space (see Fig. 2). In 
particular, when n — I, Ea. (|34p reduces to 



In Fig. 2, the behaviour of WF distributions of photon-added thermo state are plotted in {q,p) phase 
space and (|q;| ,0) space. From Fig. 2, one can see clearly the modulation action of photon-added 
number and temperature. The "oscillating frequency" of WF increases with the increasing photon- 
added number; while the amplitude value of WF in (|a| ,9) space decreases with the increasing 
temperature (corresponding to 9), which indicates that the nonclassicality is weakened at finite 
temperature. 



Uisng Eq.(I32D and the WF W2 (a) of p2 is given by 



W2{a) = C2 (0, 0| S-t (6*) a"5 (6*) [5^ (6*) A (a) S* (6*)] 5^ (61) a^"5 (6*) |0, 6) 
= {n,0\S'' {9)A{a)S{9)\n,0) 

_ (_l)"e~2|apBcch29 ^ n! f 4cosh'g , 2V 




(34) 




(35) 



VI. WIGNER FUNCTION OF THERMO NUMBER STATE 



At finite temperature, according to TFD, the number state \n) is replaced by \n,n) , thus the 
thermo number state (i.e., number states at finite temperature) is S (9) \n,n) in the enlarged Fock 
space. Using the un-normalized coherent state representation of number state. 
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where \z^z) = exp[za^ + za^ ] |0,0) is the non-normahzed two-mode coherent state, and employing 
Eg. p^ . we calculate the WF W3 (a) of thermo number state as 

Wsia) = {n,n\S^A{a)S\n,h) 

= A-—, ;— exp{-(te+ /r)tanh20 

+ (ri - /z) sech26l + zS* + + ri^* + tF}^^^^^^,^o , (37) 

where we have set 

g-2|Qpscch2e 

A = — ^ — , E = 2asech26' cosh 9, F ^ 2asech26' sinh 6'. (38) 

Trn!"' cosh 20 

Expanding the exponential term exp [{rt — fz) sech2^?] as series, we have 



X E ^-^^-^^^irtyif^)'^M^^E* + fE + tF + rFX^,^^^^^, 



(-1) sech'+''26l 

Ukl dF^dF*^ dE'^dE*'' 

Lk=0 



J2n J2n 

^ df^ d?^ [- + ''^ 20 + /i? + rF* + zE* + tFl^,^^^^^, . (39) 
Then making the variable replacement for /, r, t, z we can rewrite Ea. (j39p as 



Wsia) = y^tanh2"26l ^ 

i,fc=0 

^ dp^dr"^ dz^dt^ 



32/ 



ak 



Ukl 



dF'dF*^ dEf^dE 



^tanh2"26l ^ 



exp —fr + JE 
(-l)''sech'+^'20 



rF* 



tanh 20 



-tz + zE* 



tF 



tanh 29 



z=t=f=r=0 



l,k=0 
Qk+l Qk+l 

^ dEi'dF*^ dE*^dFi 
Noticing the formula 

Ql+k 



m 



Hn,n E, 



F* 



tanh 29 



Hn,n E 



F 



tanh 29 



Q^lQj^k 



{m-l)\{n~k)\ 



HjTL—l.n—k 



we have 



^,2g-2|a|^scch2e ^ (_!)'= sech'+^20 tauh^'"-') 20 



TT cosh 20 ^j^^ l\k\[{n-l)\{n^k)f 



Hn~k.n-l E 



p* 



tanh 20 



(40) 



(41) 



(42) 



From Eg. p^ one can see clearly that the WF of thermo number state is a real number. 
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In particular, when n = 0, noticing that tanh6' = e~2'^P^ cosh 9 = i_^-ri^ , sinh 6* = j^rj^^, 
Eq.(I32) reduces to the WF of thermo vacuum state |0(/3)) in Eg. ipD]) . On the other hand, when 
T 0,(i.e., finite temperature case reduces to zero temperature case) e~^" — > e^°° — > 0, sinh^? 0, 
cosh 61 I, E ~* 2a, ^.^^^ a*, and noticing Eg. ([29]) and the definition of two- variable Hermite 
polynomials [24,25], 



1=0 

which leads to Hn-k,o (Sa, a*) = (2a)"^'^ , then Eq. (|42)) becomes 



A:! [(n- fc)!] 

- (-i)" .-2ior— 



^„fc![(n-fc)! 

= ^e-2|"I^L„(4|ar), (44) 

TT 

which is just the WF of number state |n) at zero temperature. 

In sum, by using TFD and Weyl ordered operators' order-invariance under similar transforma- 
tions, we present a new approach to deriving the exact expressions of Wigner functions for thermo 
number state, photon subtracted and added thermo vacuum state. These WF are related to the 
Gaussian-Laguerre type functions, which are easily to be further analysed. The affection of tem- 
perature to nonclassical behaviour of the fields is manifestly shown. For discussions about the 
decohercnce at finite temperature, we refer to [30,31]. 

Appendix A Checking Eq. ipp)) 

In fact, in original Fock space, the photon-subtracted thermo state is expressed as [20] 

pi - CiTra [a" lO(/3)) (0(/?)l a^"] = Cia"pea^", (Al) 
where pc is the thermo state 

oo ; -. 

^^ = Er^TTITTlO(^l = ^e''^"'"^. nc = ^inh'e. (A2) 
1=0 [nc + 1) "'^ + ^ 

Using the the coherent state representation of Wigner operator [26] , 

/_,2 
-^lz)(-z|exp[-2(0a*-z*a)], (A3) 

where \z) is the coherent state [27,28], we have 
lyi(a) = Tr(A (a) pi) 



Note that 



and 



- / ^(-0|a"e"'"'"^a^"|z)exp[-2(za*-z*a)]. (A4) 



[ric + 1) 



(A6) 
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Eq.(A4) can be rewritten as 



Wi [a) 



(n, + 1) 
X exp 



JlcZ 



2 — 2 (za — z a) 



+ 

Further using the operator identity [29] 

aV" = (-l)":i/„,„(zat,za) : , 
where i?m,n (^^y) is the two-variable Hermite polynomials, whose generating function is 



Hm,n {x, y) 



exp [-tt' + tx + t'y] 



t=t'=0 



we have 



Wi{a) 



.exp<i-^^^^|z|' + f .11!^ _ 2a* ) z + (2a - ^^) z' 



exp 



(-1)" e^l"l' a^" 
nlTT {uc + 1)" 2nc + 1 df^dt 
ric + 1 / irn, 



Uc + 1 



■ exp [— ir] 



t = T=0 



2nc + 1 \ ric + 1 



- 2a* (2a - it) 



exp 



t=r=0 



n!7r(nc + l) 2nc + l9t"5T" t 2nc + 1 
2ia* ^ t + 2ia 



-tT 



ric 



2nc + 1 2nc + 1 

e~OT (-If ^ 



n+l I Hn.n 



t=r=0 

2mra 



(2nc + l) ^ 
then using the relation 

(-1)" 



,2i 



y(2n, + l)(n, + l)' V2n, + 1 



ric + 1 



-a 



i7„,„ (x,y) = L„ {xy) , 



and noticing that ric = sinh^ 0, 2nc + 1 = cosh 26*, Eq.(AlO) can be put into 



Wi (a) 



e / 4?!^, |a|' 



7r(2ne + l)"+' " V 2nc + iy' 
which is just the Ea. ([5(I| . Thus we have checked the result using a new appraoch. 



(A7) 



(A8) 



(A9) 



(AlO) 



(All) 



(A12) 
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